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We consider evolution equations of the Lotka-Volterra type, and elucidate especially their formu- 
lation as canonical Hamiltonian systems. The general conditions under which these equations 
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related to the Liouville model on a lattice, is considered in detail, both as a classical and as a 
quantal system. 
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1 Introduction 



The Lotka-Volterra model fT], H is defined by a set of non-linear differential equations of the 
following form, 

7 N 

= eAWA+J2 FabWaWb , A= 1,2,.. .N (1) 

<^'' B=l 

In the equations above the quantities wa, ^ = 1,2, ...,N are to be determined in terms 
of the given constant quantities and Fab- The Lotka-Volterra equations are principally 
known as qualitative models in population dynamics with interacting species, as well as rate 
equations for chemical reactions with constituents. However, the equations in question have 
actually turned out to be rather universal with different applications in physics , as witnessed 
e.g. by the case of the Liouville model formulated on a lattice 0. 

The matrix F in Eq.(0) ought to fulfill certain conditions which relate to the concept 
"crowding inhibits growth". This condition will here be related to the requirement that the 
matrix F be antisymmetric. 

Fab = —Fba (2) 

We wish to elucidate the circumstances under which the equations including the anti- 
symmetry condition (H), admit a canonical Hamiltonian formulation. This question has been 
answered in a satisfactory manner by Kerner already quite some time ago, however with 
the restriction to an even number of species and with the condition that the matrix F be 
regular. Here we give a general analysis, valid for even or odd A^, as well as for regular or 
singular matrices F. In the singular case (which always occurs if A^ is odd), the equations (|l]) 
admit additional conserved quantities besides the Hamiltonian, under certain circumstances. 
We analyse these circumstances in detail. In special cases this phenomenon has been noted in 
the literature (e.g. [^), and has been referred to as a multi-Hamiltonian structure. 

The analysis of the general Lotka-Volterra equations (|l]) including the condition @ is given 
in Sec. II below; the next Sec. HI then deals with a particular example corresponding to an 
odd number A^ = 3 of species. 

In Sec. IV we consider a specialisation of the Lotka-Volterra model, which appeared in 
a formulation of the Liouville model on the lattice IQ, the quantal formulation of which was 
further considered by Volkov 0] who introduced the name "Quantum Volterra model" in this 
connection. 

The Hamiltonian proposed in |^ and in ^ is actually one of the additional conserved 
quantities that appear as a consequence of the special features of the model under consideration, 
and differs from the Hamiltonian derived in this paper. Our alternative formulation permits a 
simple canonical quantization of the model under consideration. 

The final Sec. V contains a summary and discussion. 
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2 The General Lotka-Volterra Equation 



In analysing Eq. (|l|) it is convenient to introduce new variables as follows, 

= logWyl 

Then Eq. (|T|) can be rewritten as follows, 



N 



iA = (^A+J2 ^AB exp 



(3) 



(4) 



B=l 



In order to analyse Eq. (|ID further, it is convenient to refer to the socalled normal form 
of the antisymmetric matrix F. It is well known 0] that, by making an appropriate basis 
transformation, any antisymmetric N x N matrix F can be transformed into the following 
normal form: 



/ -h 
+ki 



V 



...\ 





-\-kn 



kri 





(5) 



J 



where the (positive) quantities ka, a = 1, ...,n are the square roots of those characteristic values 
of the matrix — F^, which are different from zero. The number n of non-zero characteristic values 
/c^ is thus given by the rank {2n, say) of the matrix F. 

It is furthermore convenient to relate the normal form (|^) to the following linear equations, 

N 

^ FabXcB = + kaVaA, Oi=l,...,n (6) 



B=l 



N 



FABVaB = - kaXaA, a = 1, .. 



,n 



(7) 



and 



B=l 



N 



E F^ABZ^B = 0, /3 = 1,...,N -2n 



B=l 



with the understanding that Eq. @ is empty if the rank of the matrix F is N (so that N = 2n) , 
in which case the matrix F is regular, 



det F =^0. 



(9) 



The cases even or odd differ qualitatively in general, since in the latter case the matrix 
F is necessarily singular. Then Eq. (^ has an odd number of nontrivial solutions. Eq. (H) 



2 



may of course also have non-trivial solutions if N is even, in which case there is necessarily an 
even number of such solutions. 

We ortho- normalize the vectors XaA, VaA, and Z/^a, properly, 

{Xa, Xfs) = {Ua, Vp) = {Za, Zp) = S^p (10) 

where the inner product {u, v) of any vectorlike quantities ua and va is defined as follows, 

TV 

{u,v) = '^aVa (11) 
a=i 

The equations (||), (^ and (§) imply the following orthogonality relations, 

{xo^^yp) = 0, = 1, ...,n; (12) 

and 

(x„, zp) = (na, zp) = 0, a = 1, ...,n,(3 = 1, - 2n (13) 

After these preliminaries we return to the equations (|I]). Contracting Eq. (|I]) with any 
solution Zf3 of Eq. (P), one obtains, in view of the antisymmetry condition (Q), 

{Z(,,0 = {Zf3,e)=rp (14) 

where we use the inner product notation given above in Eq. ([TT|). Thus, 

{Zf3, 0=rpt + Kp, P = l,...,N-2n (15) 

where the quantities K/3 are arbitrary constants of integration. The equations (0) define N—2n 
constraints among the variables ^a] the number of these constraints is equal to the number of 
linearly independent solutions (if any) to the zero-eigenvalue equations (^. In particular, if 

rf, = {zp,e)^0/3 = l,...,N-2n (16) 

then the N — 2n quantities {z/3,^) are linear functions of t. 

The equations (|ip still contain 2n unknowns; one obtains a convenient set of 2n equations for 
these unknowns by contracting the equations (|l|) with the solutions Xa and ya of the equations 
(H) and (^, respectively. Thus, 

N 

{Xa, i) = iXa,e) - J2 KyaB exp (17) 
B=l 

and 

TV 

iVa, i) = (?/a, e) + ^ k^x^B exp iB (18) 

B=l 
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It should be noted that the equations (|T^, and (|1^) are equivalent to the original 
equations (|I]) under the condition since the ortho normal vectors Xa, Va and Za form a 
complete set in dimensions. Thus, 

n N-2n 

U=J2 [(^"' Ox»A + iVa, OVca] + J2 O^PA (19) 

a=l p=l 



The even number of equations ( ]T7| ) and are in fact Hamiltonian equations in coordinates 
Qa and momenta pa defined as follows, 



Pa — Cai^Xa, (^) 



and 



(20) 
(21) 



where the quantities Cq, and are constants, which satisfy the following condition, 

CaD^ka = 1 (22) 

Namely, let 

N n 



H{p, q;t)^ J2 {Cb{P, q;t)} -J2 [Ca{Xa, e)qa - DaiVa, ()Po 



(23) 



B=l 



a=l 



where the quantity ^(p, q; t) is expressed in terms of the quantities Pa, qa defined above, as well 
as in terms of the quantities and Kj^ defined in Eq. (0) and (^) , respectively, according to 
Eq. (0), 



N-2n 



U{P, q;t) = J2 [i^a ^Pa)XaA + {D^^qa)yaA\ + J2 i^f^^ + Kfl)Zf3A 



(24) 



a=l 



Using Eq. (p^ and Eq. (22) one obtains straightforwardly the following results from Eq. (^31), 



and 



OH 

dpa 

dH 

dqa 



N 



J2 f^a^»B exp {Csip, q;t)} + {ya, e) 

.B=l 



N 



J2 ^aVo^B exp {^b(p, q]t)} - (a;^, e) 

.B=l 



(25) 



(26) 



Then, using the relations just given, it is a simple matter to verify that Eq. ( p!7D is nothing but 
the following. 

Likewise, Eq. (|18|) is equivalent to the following. 



dqo 



+ 



dH 



dt ' dpa 
4 



(2^ 



It should be noted that the Hamiltonian (^) is in general explicitly time- dependent if the 
zero-eigenvalue equations have non-trivial solutions. This requires that the matrix F ocur- 
ring in the Lotka-Volterra equations (P be singular, which is always the case if the dimension- 
ality N of the system is an odd integer. Be that as it may, as we have just demonstrated, there 
is in any case always an even- dimensional subset of the equations (|l]) which can be written in 
canonical Hamiltonian form, in terms of canonical momenta (pO|) and coordinates (|2T|), with 
the expression (^) as a Hamiltonian. 

The explicit time-dependence of the Hamiltonian (^) disappears if the vector e (the array of 
rate-constants ca) occuring in Eq. (|I]) is orthogonal to all the solutions of the zero-eigenvalue 
equation (H), 

{Zf3,€) = , P = l,...,N-2n (29) 
If the equations (|29| ) are in force, we have, according to Eq. (|15|), 

{zf,,0 = Kp, f3 = l,...,N-2n (30) 

Each orthogonality condition (p9|) with a fixed value of the index P, gives rise to a conserved 
quantity, namely the corresponding expression in Eq. (0), according to Eqns. ([I5|). 

We conclude this section by summarizing the results obtained so far: 

The Lotka-Volterra equations (||), including the condition (0), can always be written as 
a system of canonical Hamiltonian equations in an even number of unconstrained canonical 
variables, together with a set of explicitly solvable constraints, which are time-dependent in 
general. This set is non-empty, i.e. the constraints in question occur in general, if the matrix 
F in Eq. (|l]) is singular. The number of constraints equals the number of linearly independent 
eigen-vectors z/s of the matrix F, which correspond to zero-eigenvalues. 

For non-singular matrices F the Hamiltonian does not depend on time explicitly. In the sin- 
gular case, the Hamiltonian is in general explicitly time-dependent. This time-dependence dis- 
appears from the Hamiltonian, if all the eigen-vectors zp are orthogonal to the vector (ei, eAr) 
made up by the rate-constants in the basic equations (|ip. Each orthogonality condition of the 
aforementioned kind gives rise to a (time-independent) linear constraint among the variables 

in the equations (^), which are equivalent to the basic equations (|T[). 
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3 An Example Involving Three Degrees of Freedom 



We consider below a system of the Lotka-Volterra type, with three degrees of freedom, which 
has been analysed in great detail by Grammaticos et.al with the objective of making a 
systematic search for first integrals of the system in question. A special case of this system has 
been presented by Nutku |0 as an example of a bi-Hamiltonian system of the Lotka-Volterra 
type. 

The general system is the foUwing (a, b, c and A, fi, v are constants), 

s = s(A + ct + u) ,t = t{fj, + s + au) , u = u(z/ + bs + t) (31) 

Rescaling the variables in Eq. (|3TD as follows, 

111 
s — > — s = wi ,t — t = W2 u —' — u = W3 (32) 

ai OLi Q!3 

where the parameters a„,r2, = 1,2,3, are at our disposal, and using the ^-variables introduced 
in (H), we get an equation of the form (H) from Eq. (pTj), 

3 

= + XI ^^-B exp is (33) 

B=\ 

where, 

ei = A , £2 = /i , ez = v (34) 

and 

/ a^c as \ 
F = a\ a^a (35) 
\a\b a2 J 

The matrix F defined by Eq. (|35D does not fulfill the antisymmetry condition (^ as such, but 
becomes antisymmetric under the foUwing conditions, 

CKi = —a2C , a2 = —a^a , 0:3 = —aib (36) 

But Eqs. (^) can be fulfilled if and only if 

abc = -1 (37) 



From now on, we assume that the parameters a, b, c satisfy the condition (|37|). Without essential 
loss of generality, we can choose one parameter a„ at will in the equations above. Taking ai = 1, 
we satisfy the equations ( p6| ) as follows, 

ai = 1, a2 = ab, = — b (38) 



We then apply the results of Sec. II to the system of equations above. 
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In the first place we have to consider the following system of linear (eigenvalue) equations 
(compare (||), (0) and (g)), 

3 3 3 

FabXb = + ki/A , PabVb = - kxA , Y ^abZb = (39) 

B=l B=1 B=1 

The following are appropriately orthonormalised solutions to these equations, 



? ^ if- - 1 1 



/ ab 

-b I (40) 
V 1 



with 



k = Vl + &2 + a262 (41) 
Using Eqns. (^) one now obtains the follwing constraint, according to Eq. (|15|), 

{z,0 = l [ab^i - hi2 + 6] = rit + K, (42) 



n = (2, e) = ^[a6A-6/i + z/] (43) 



where 



The remaining variables can then be taken to be the canonical variables defined in Eqns. ([20| , 
|2ll) , i.e. in the present case, 

p^C{x,i) = ^£=={i, + bi,) (44) 

and 

q ^ D{y,0 = - , /^^ ((l + + ^^''^ - aft^a) (45) 
where the constants C and D are related as follows 

CDk = 1 (46) 

but are otherwise arbitrary. 

The equations (^2]) , (^if) and (|45|) can easily be inverted. 



6 = ^(rit + irO-v/TTfe^Cg 

1 / ^ bk cib ^ 
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Using the results above, it is simple to verify that the variables p and q defined in the 
equations (|4^), (|45| ) above, satisfy canonical Hamiltonian equations, i.e. the equations (pTj), 



(psD , with the following Hamiltonian i/, 



1 3 
H=—[{^i + hv)ii - (A - ahv)£,2 - h{\ + a/i)e3] + E ^xp (48) 



yl=l 



The Hamiltonian (^) is in general explicitly time dependent, since the sum of exponentials 



in (^Sf) above, depends on time in general, as can be inferred from the equations (|4^) . However 



if the orthogonality condition p9| ) is in force, i.e. if 



ri = (z,e) = i[a6A-6// + z/] =0 (49) 



then the Hamiltonian ( ^HD becomes a constant of motion. However, if the condition p9|) is in 



force, then there exists a second constant of motion, according to Eq.(|4^), namely the following, 

\[ahi,-bi2 + i,] (50) 

The appearance of the two conserved quantities (^) and (0) under the condition ( ^Uj) is 
just the phenomenon which has been referred to by the term bi-Hamiltonian in the present 
context by Nutku H. 



This terminology is perhaps a little unfortunate, despite the fact that the quantity (|50D can 



be understood as a Hamiltonian for the system (|33[) (if condition (|9) is in force) in a certain 



general sense ||T0[- On the contrary, the Hamiltonian (^q) , is a Hamiltonian in a strict sense, 
regardless of whether the condition ( ^9]) is true or not. By strict sense is here meant that the 
formulation involving the Hamiltonian (^) is explicitly canonical, involving a known pair of 
canonical variables p and q. So, in the general case, only the quantity (^8]) merits a designation 
as a Hamiltonian in the strict canonical sense employed in this paper. 

The considerations above are naturally generalisable to systems with more than three de- 
grees of freedom, in view of the general results given in Sec. H. 
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4 The Liouville Model on a Lattice 



The Liouville model on a lattice discussed by Faddeev and Takhtajan and by Volkov 
is specified by the following classical equations of motion, which is a special system of Lotka- 
Volterra equations, 

dwA I , 

—wa{wa+i - wa-1) , A = 1,...,N (51) 



dt 2A 

where is an even positive integer, and where A is a parameter. The variables wa are 
furthermore supposed to satisfy the following periodicity condition. 



WN+n = Wn , n = 0,1,2, 



(52) 



Comparing Eqs. (|5l|) with the general Lotka-Volterra equations (|l]) one observes that the model 
defined above by the Eqns. (^) corresponds to the case of vanishing rate constants eA, 



eA = 0, A=1,2,...,N 
and to the case in which the matrix F has the following form. 



(53) 



/ A 

-A 



(Fab) 



V A 



A 



-A 



A 
-A 



A 
-A / 



(54) 



where 



The matrix (Q) is singular, 



det F = 



(55) 



(56) 



This circumstance together with the conditions ( ^31) implies the existence of additional con- 
served quantities (constants of motion) besides the Hamiltonian, according to the general results 
presented in Sec. II. 

At first we give a brief analysis of the canonical formulation for the classical latticised 
Liouville model, along the lines presented in Sec. II. 
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4.1 Canonical Hamiltonian Formulation of the 
Classical Liouville Model on a Lattice 



First we rewrite Eq. (|5lD using the variable (|^) 



N 



= A (exp ^^+1 - exp = ^ Fab exp (57) 



B=l 

where the matrix F given by the expression (^^. According to the general method developed 
in Sec. II, we consider to begin with the eigenvalue-equation (||), i.e. the following equation, 

N 

E FabZ(3,b = (58) 

B=l 

As is immediately verified, Eq. ( pHD has two linearly independent (normalized) solutions Zf^, 
which we take to be the following, 

^i,A = ^(1 - (-1)^) , A=1,...,N (59) 

and 

^2,A = -^(1 + (-1)^) , A=1,...,N (60) 



The existence of the two non-trivial solutions (|59D and (|60D implies the existence of the 
following two conserved quantities, according to Eqs. (|^) and (|30D , 



and 



|E?2A-l = i^l (61) 

^ i: 6a = K, (62) 

A=l 

Needless to say, the Hamiltonian (pSf), which now takes the following simple form, is also 
conserved, 

N N 

H = exp ^B = wb (63) 

B=l B=l 

We then have to consider the equations (^ and (0) with the matrix F given by (|5^ . 
After calculations one finds the follwing solutions XaA and i/aA, 

x^A = ^i=(l + (-1)^) [cos(A(^„) + sin(A^„)] , a = 1, 2, - 1. (64) 
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and 



where 



VaA = -^(1 - (-1)^) [cos(yl^„) - sin(A(^„)] , a = 1, 2, - 1. 



= , a = 1,2, ...,-iV - 1. 



iV ' ' ' '2 

The eigenvalue parameter is expressed in terms of ip^ as follows, 

ka = 2A sin (^9^ , a = 1, 2, -A^ — 1. 



(65) 
(66) 

(67) 



Applying the general formulae ( PP] ) and (^Tj) we now obtain the following expressions for 
the canonical momenta and coordinates, respectively. 



Pa 



and 



^;J^i:6A-icos[(2A-i)^„+- 



(68) 



(69) 



A=l 



The inverses of Eqns. (^Sj), (^) are obtained straightforwardly, either by using the general 
result (^) or by evaluating the appropriate sums involving trigonometric functions. 



a a 13=1 



(70) 



where the constants are given in Eqns. (|6TD , (|62D. 



We conclude this sub-section by giving a Poisson-bracket formulation of the equations (pTl). 
The basic Poisson-bracket is the following 



{^a,^b}pb = 



a=l 



dqa dpa dqa dpa 



(71) 



Finally, evaluating the relevant trigonometric sums in ([7T| ) (compare Eqns. (p^), (|65|)), one 
obtains, 

{^a,^b}pb = ^{^A,B-i — Sa,b+i) = Fab (72) 



Using the result (^) it is straightforward to verify that the equations (|57|) are reproduced 
by the following expressions. 



N 



iA = {Ca,H}p^ = J2 {Ca,^c}pb 



dH 



C=l 



c 



(73) 
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where the Hamiltonian is given by Eq. (|& 

Reverting to the original variables Wa (compare Eq. (^), one gets the following Poisson- 
bracket relations from (ff^), 

{wa, wb}pb = XwaWb {Sa,b-i - Sa,b+i) (74) 



It is perhaps necessary to emphasize that the result (ff^) is firmly based on the canonical 
structure derived above for the latticised Liouville model, and not merely an example of a 
bracket structure which is consistent with the Hamiltonian (^Sj) and the equations of motion 

In the aforementioned paper by Volkov [0, an alternative Hamiltonian formalism to the one 
presented above has been proposed for the quantization of the system defined by the equations 



(HD- 

The Hamiltonian Hy used by Volkov is simply the sum of the two constants of motion (pTl) 
and (|62D, apart from a multiplicative constant, 

1 ^ 

A=l 

where 7 is an arbitrary coupling constant, which is introduced for convenience. The quantity 
(|75D is a constant of motion, since it is the sum of two constants of motion (compare Eqns. (|6l| . 



The Hamiltonian scheme of Volkov consists of a bracket formulation of the equations (|51|) 
with (^) Hamiltonian. Thus, 

wa = {wa,Hv} = y2 {'^a,wb} 1:: (76) 

The basic bracket given by Volkov is the following (after correction of a crucial printing error; 
an overall factor waWb is missing from Volkovs expression, Eq. (23) of Ref. [0), 

1 

{wa, wb] = -IWaWb [(4 - Wa- wb){Sa+i,b - Sa-i,b) + wa-iSa-2,b - wa+iSa+2,b] (77) 



It is indeed simple to verify, that the equations (|51|) are reproduced as bracket equations 
(|76D if one uses the bracket ([77|). However, it should be noted that the Hamiltonian Hy defined 



in Eq. ( [75D which is used in this context, is not bounded from below, and does therefore not 
permit the construction of classical statistical mechanics for the (classical) system in question. 
On the other hand, the Hamiltonian H defined in Eq. (0), which has been obtained here on 
the basis of general considerations of the canonical structure of systems of the Lotka-Volterra 
type, is bounded from below, and is thus preferable also from this point of view. 
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4.2 The Canonical Quantization of the 
Liouville Model on a Lattice 



The quantization of the model considered above is straightforward if one uses the canonical 
structure derived here as a starting point. Following the usual prescription of replacing Poisson 
brackets involving classical canonical variables by commutators involving the corresponding 
operators one has the following canonical quantum operator relations, 

[Pa^Pfi] = 0, [qa^q/s] = 0, [qa,Pl3] = ih6a(3 (78) 

From the canonical commutators ( [78l) follows indeed the quantization prescription given below 
in terms of the basic variables ^ of the problem, 

{U,^B}pB — ^[iA,iB] (79) 

where quantal variables are distinguished by a hat, thus: C,. 

The Hamiltonian H given in Eq. (|63|) can immediately be generalized to an operator Hamil- 
tonian H. Using the inverse relations ([fOD, or more explicitly the following relations. 



4A ^"^"^ TT fJ 

6a = -7^ E sin sm[{2A)ipa + j\ + V T?^^ (§0) 

V ' a=l ' 



1) 



and 

4A TT / Z 

6a-i = E '^"^^ sin cos[(2yl - l)^?^ + 7] + V T7-^i 

ViV Q,=i 4 V -'V 

one obtains the following expression, 

H = Tip) + Viq) (82) 



where 



and 



N 
2 



2 4 A ^ "'^ TT 

T{p) = E exp <! \lj^K2 + -7^ E D^p^smip^sm[{2A)ip^ + ^] \ (83) 



A=l » ^ ' V iv 



JV 

2 



Y 4A ^^"^ 

y(g) = Eexp<|^/^i^i + -^ E agaSin</.,cos[(2A- !)</.„ + J ^ (84) 



A=l » ^ ' V iv 



The Hamiltonian (^) is expressed as a sum of terms depending separately on the momentum 
variables p and coordinate variables q, respectively, so there are no essential quantum ordering 
problems in the quantization procedure outlined above. 

The quantization procedure considered by Volkov ^ is not explicitly canonical and can 
therefore not be compared directly with the canonical standard procedure given above. 
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5 Summary and Discussion 



In this paper we have given a general analysis of the canonical structure of the system of 
differential equations which are known as the Lotka-Volterra model, supplemented with an 
antisymmetry condition (which is crucial for our analysis) which is also frequently assumed 
in connection with the Lotka-Volterra equations. Altogether this defines a dynamical system, 
with a finite number of degrees of freedom, which we call the Lotka-Volterra system. 

It has been shown that the Lotka-Volterra system always can be resolveded into an explicitly 
canonical Hamiltonian subsytem, involving an even number of canonical equations, together 
with a set of explicitly solvable constraints, which are time- dependent in general. The set of 
constraints is empty if the matrix which defines the interaction between the various components 
in the model, is regular. A set of canonical variables (pairs of coordinates and momenta) is 
explicitly constructed for the Lotka-Volterra system. 

Furthermore, it has been shown that the Lotka-Volterra model gives rise to several conserved 
quantities (constants of motion) if the parameters of the model satisfy certain conditions, which 
have a general geometric characterisation as certain orthogonality conditions. The occurrence 
of such additional constants of motion has been termed a "multi-Hamiltonian structure". 

As an illustration, a Lotka-Volterra system involving three degrees of freedom has been 
considered in detail. 

An important special Lotka-Volterra system, which is associated with the Liouville model 
on a lattice has also been analysed in detail. The canonical structure and Hamiltonian formu- 
lation of the system in question has been shown to be a straightforward consequence of the 
general formalism developed in this paper. The canonical quantization for this case, which 
follows straightforwardly from the canonical structure of the underlying classical model, has 
been outlined in detail. This canonical quantization procedure differs in a non-trivial way from 
a previously proposed Hamiltonian procedure |^ for the model under consideration. 
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